Stylized Facts Generated Through Cellular 
Automata Models. Case of Study: The Game 

of Life 



H.F. Coronel-Brizio, A.R. Hernandez-Montoya *, 
M.E. Rodriguez- Achach and G.A. Stevens-Ramirez 

Facultad de Fi'sica e Inteligencia Artificial. Universidad Veracruzana, Apdo. 
Postal ^15. Xalapa, Veracruz. Mexico 



Abstract 

In the present work, a geometrical method to generate a two dimensional random 
walk by means of a bidimensional Cellular Automaton is presented. We illustrate 
it by means of Conway's Game of Life with periodical borders, with a large lattice 
of 3000 X 3000 cells. The obtained random walk is of character anomalous, and its 
projection to a one dimensional random walk is analyzed, showing that it presents 
some statistical properties similar to the so-called stylized facts observed in financial 
time series. We consider that the procedure presented here is important not only 
because of its simplicity, but also because it could help us to understand and shed 
light on the stylized facts formation mechanism. 
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1 Introduction 

Recent methodologies from Physics have been successfully applied to the study 
of Economy and Financial Markets Complexity. This new area of research is 
called Econophysics. The econophysics community is currently doing active 
research in topics such as the study of the distributional properties of the 
variations of the Stock Market, Network Analysis of Economical Phenomena, 
Financial Crashes, Wealth Distributions, study of Financial Markets under a 
Microscopic point of view, etc. About the latter, it involves techniques named 
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Microscopic Simulation (MS) [1] , that consists of the study of a complex sys- 
tem by individually following on a computer each agent of the system and 
its interactions with the other agents, and simulating the overall system evo- 
lution. This technique is becoming very promising and useful. Following this 
direction, very interesting models that reproduce the statistical properties of 
financial Markets, named by the economics community "Stylized Facts" [2,3], 
have been proposed [4,5]. Also in a very related way, a great amount of work 
has been carried out by constructing artificial stock markets by means of Cel- 
lular Automata models [6,7,8]. 

Cellular Automata (CA) are spacetime-hke discrete deterministic dynamical 
systems whose behavior is defined completely in terms of local interactions. 
Cellular automata were introduced by John Von Neumann who intended to 
understand the biological mechanisms of self-reproduction [9] and have now 
become an important object of study because of their intrinsic mathematical 
interest as well as their success as a tool to model complex phenomena in 
physical, chemical, economical and biological systems, design of parallel com- 
puting architectures, traffic models, programming environments, etc [10]. 
Cellular automata became very popular at the beginning of the 70's thanks 
to an article written by Martin Gardner and published in Scientific American 
[11], about the cellular automaton called The Game of Life (GOL) or just 
Life. This cellular automaton was invented by the mathematician John H. 
Conway at the end of 60's and since, it displayed a very rich and interesting 
complex behavior; very soon it became the favorite game of the -at that time 
incipient- community of computer fans. 

LI The Game of Life 

Life is a class IV bi-dimensional (shows complex behavior) totalitarian cellular 
automaton [12]. The updating rule that determines the Game of Life evolution 
is applied on a Moore neighborhood as follows: a) a dead cell surrounded by 
exactly three living cells is born again, b) a live cell will die if either it has 
less than two or more than three living neighbors. This simple rule produces a 
very rich and complex behavior, generating self-organized structures and also 
producing very important and interesting emergent properties (formation of 
self- replicating structures. Universal Computation, etc). It is for this reason 
that we have chosen Life to generate a bi-dimensional random walk and then a 
time scries with expected complex characteristics. In section 2 we explain how, 
by using a geometrical procedure applied to Life, we can generate a time series 
Tj, analyzed in this work, that shows statistical properties very similar to those 
of Financial Time Series, known in the economic community as stylized facts. 
In section 3 we show the numerical results corresponding to these analyses. 

2 Generating a random walk by the Game of Life. Defining our 

observable 

Due to the need of dealing with a finite size N x N of the CA lattice, in our 
work we have chosen to program our CA using periodical boundary conditions. 
After centering a Cartesian coordinate system XY in the lattice of the CA as 
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shown in figure la), in order to obtain the one-dimensional observable analyzed 
in this paper, we construct the position vector R{i)cM of a point following a 
bidimensional random wallQ obtained for each time step of the M generated 
points, i = 1,2,3,.., M, as follows: 

^ ' 1 > 5 5 5 ^ N N 

R{1)CM := {Xit)cM,Y{t)cM) = EE (1) 

k=l 1=1 

where N x N,is the total number of cells of the CA, Cki{i) denotes the state (1 
or 0) of the cell in the coordinates {xk{i),yi{i)) with respect to the coordinate 
system XY at time step i. So, while the CA evolves, this vectoiO] is calculated 
for each time step i = 1,2, 3, .., M. Figure lb) shows 10000 time steps of the 
evolution of the vector R{i)cm- Our observable will be the distance to the 
origin of R{i)cM, i-e. := ^X{i)lj^j + Y{i)l;M, i = 1,2,3,..,M. In order 
to analyze the one-dimensional time series generated by as time goes on, 
we construct the returns or logarithmic differences Si, of the time series rj as 
usual: Si := logr^+i — logr^, for i = 1,2,3,..,M. 




Fig. 1. Left panel: Coordinate system used to calculate the vector R{i)cM for each 
state of the GOL CA (only for illustration purposes). Right panel: Random Walk 
generated when R{i)cM evolves in time. 



2.1 Data Sample 

In order to have certain control of GOL initial states, configurations of 20%, 
40% and 60% of living cells were initially set up at random on lattices of size 
3000 X 3000. For each initial density, 20 random walks of 20000 steps were 
generated, giving us a total 1200000 steps to analyze. Considering certain 
characteristics such as the finite size of the GOL lattice, or a particular initial 
configuration, GOL generated fiuctuations tend to die out after an undeter- 
mined number of time steps. Because of this we have applied a cut off to r^, 
i = 1,2,3, ..,M, considering its values such as M < 3000. Figures of next 
section will show the reason of this data selection more clearly. 

^ This random walk is of character anomalous, but this is demonstrated in a forth- 
coming paper. 

^ Any resemblance with the Center of Mass (CM) definition of a mass distribution is 
true, R{i)cM describes how the "center of mass" of living cells distribution evolves. 
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3 Numerical Results: Statistical Properties of returns 

In this section we will show that the time series rj, z=l,2,3,..,M has statis- 
tical properties that resemble the stylized facts of financial time series. Figure 
2 shows a plot of for three different time ranges, while figure 3 shows the 
graphic of returns vs. time for two different time ranges. Both figures were 
obtained from one of the typical generated random walks. 




Fig. 2. Time evolution of our observable, R{i)cM distance to Origin for time ranges 
of a) 20000 b) 8000 and c) 400 time steps. 




Fig. 3. Returns evolution. In a) is clear the effect of the CA fluctuations dying out. 
b) Zooming in the upper figure. We can appreciate the clustering volatility effect. 



3.1 Returns Distributions 

Returns distribution Si, for all of our generated samples are shown in figure 
4a). Figures 4b) and 4c) show the returns distribution in a log- log plot for 
right and left tails respectively. After applying the cut off to the sample, we 
construct the Probability Distribution Function (PDF) of the returns distri- 
bution; this is shown in figure 5. As it can be seen, both PDF tails decay as a 
power law with exponent a ~ 6; Although it seems that there are more than 
one power law present along both curves. 
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Fig. 4. Returns Distribution, a) Full Sample, b) Returns Distribution Log-Log Plot 
for right tail; c) Returns Distribution Log-Log Plot for left tail. 
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Fig. 5. Returns PDF. a) Right tail b) Left tail. In total, after the cut off we have 
analyzed 479940 events. Red lines are only sight guides. 

3.2 Returns and Absolute Returns Autocorrelations 

Figures 6 and 7 show the Auto Correlation Function (ACF) of returns and 
absolute returns for a few realizations of our experiments respectively. It can 
be seen that returns ACF shows no memory, decaying to noise level almost 
immediately (in fact figure 6 looks similar to the ACF of a daily financial ACF 
time series). On the other hand, ACF of absolute returns decays quickly to 
a positive level of about 0.2 and stays fluctuating there, showing a very long 
range memory; all this in agreement with the stylized facts. 
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Fig. 6. Returns ACF calculated for 8000 time steps and different time ranges. 
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Fig. 7. Absolute Returns ACF also calculated for 8000 time steps. Different time 
ranges are shown 

It is important to remark that although curve of figure 6 does not show corre- 
lations, a periodic pattern in each ACF there exists (although it was not shown 
here); also we have to mention that figure 7 should look different wether the 
Absolute Returns ACF were calculated along the full 20000 time steps, show- 
ing a slower decay more similar to that of financial time series. This behaviour 
was not showed here also. 



3.3 Volatility Distribution 

In this section, the volatility V{t) is calculated [13] by averaging the absolute 
returns over a time window T = nAt as follows: 

-I t+n-l 

V{t):=- E 1^(01, (2) 
^ t'=t 



Here we have set up At = 1 time lag and a window of 50 time steps. Figure 
8a) displays the volatility for the first 3000 time steps of our observable for 
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one of our generated random walks. Lower frame of same figure displays the 
corresponding volatility distribution. Figure 9 shows that volatility PDF de- 
cays asymptotically as a power law with a ~ 13, again in concordance with 
the stylized facts. Although also in this figure, as in the PDF of retuns case, 
more than a single power law are present along the curve. 




"o 500 1000 1500 2000 2500 3000 Time 




Fig. 8. a) Volatility for a typical generated random walk with a time window of 50 
time steps, b) Its distribution. 




Fig. 9. Volatility Probability Distribution Function for all the sample. It decays as 
a power law with exponent a ~ 9. 

4 Discussion and Future Work 

The objective of this work is to generate in the simplest possible way a time 
series that displays the stylized facts that are well known to the economic 
community. By means of a very simple geometrical method, and using a bidi- 
mensional Cellular Automaton, in our case the Game of Life, we generate a 
synthetic time series with the properties of having returns distribution with fat 
tails, clustering volatility, practically no returns autocorrelations, long mem- 
ory in ACF of absolute returns and volatility distribution decaying as a power 
law. We believe that the use of our scheme to generate artificial financial time 
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series can be of help in understanding the underlying mechanisms that govern 
the formation of the stylized facts arising in real markets. Finally, and about 
the future work that has to be made, is important to understand if there exists 
a multipower or a mix up of power laws in returns and volatility distributions 
as well as the patterns presented in the ACFs of the returns and absolute 
returns. 
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